gitale 


imulatie van de 


naloge 


ekenwijze 


Dit verslag is in twee gedeelten te splitsen. De delen 
I t/m VII behandelen een ontwikkelae eenvoudige inte- 
grator en enkele daarmee bereikte resultaten. 

Met behulp van het daarmee verkregen inzicht is het 
uiteindelijke DISAR-sybteem (DIgitale Simulatie van de 
Analoge Rekenwijze) ontwikkeld, zoals beschreven is 

in het tweede gedeelte (VIIl t/m XII). De nadruk ligt 
in dit gedeelte op de verkregen resultaten (XI) en op 
de eigenschappen (X) van het DISAR-systeem. 

Hoewel niet gepretendeerd kan worden dat iedere analoge 
opstelling met het DISaR-systtem gesimuleerà kan worden, 
hebben de resultaten aangetoond dat digitale simulatie 
in sommige gevallen voordelen kan bieden boven de ana- 
loge rekenwijzee 
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INLEIDING. 


De analoge en digitale rekenwijze hebben ieder hun 
eigen specifieke voor- en nadelen wanneer ze gebruikt 
worden voor het oplossen van een bepaald probleem. 


Door digitale simulatie van de analoge rekenwijze 

is een ander vatroon van voor- en nadelen te verkrijg- 

en, hetgeen in sommige gevallen de oplossing ten 

goede komt. De digitaal gesimuleerde rekenwijze heeft 

als voordelen boven de echte analoge rekenwijze: 

1, dat vermenigvuldigers en delers principieel een- 
voudiger te verwezenlijken zijn, 

2, dat het aantal elementen of de verdeling van de 

aantallen verschillende eiementen slechts door de 

geheugenruimte van de digitale rekenmachine wordt 

beperkt, 

dat logische Ee en besturing aan ae hand 


van beslissingskriteria eenvoudig te verwezenlijken 
zijn en 

4, dat de uitkomst in veel meer decimalen te ver- 
krijgen is dan in het analoge geval, (zie echter 5). 

Bezwaren van de digitale rekenwijze die meespreken 

zijns 

2e De convergentie van ae oplossingsmethode is niet 
altijd zodanig dat de verkregen cijfernauwkeurig= 

+ heid in de uitkomst significant is. 

Oe Het kontakt tussen mens en rekenmachine is niet 
zo direkt als bij de analoge rekenwijze. 

7e De rekentijd neemt toe als het aantal elementen 
toeneemt. 

De meeste aniere voordelen van de analoge rekenwijze, 

bijvoorbeeld het eenvoudig kunnen opstellen en doorzien 

van het schema van het probleem, blijven. 

Het doel van het in dit verslag beschreven werk is 

nu als volgt te formuleren: 

Het digitaal simuleren van de analoge rekenwijze, met 

inachtneming van de volgende eisen: 

a. Het mag niet nodig zijn de differentiaalvergelijking 
te kennen. Het analoge schema alleen moet voldoende 
zijn. 


Db. De rekenwijze moet zo nauw mogelijk aansluiten bij 


de analoge rekenwijze, zodat grootheden, die op 
de analoge rekenmachine als uitgang te verlangen 
zijn ook digitaal als uitgangsgrootheden kunnen 
optreden. 
C. De nauwkeurigheid moet zo mogelijk groter zijn 

dan bij analoge rekenmachines optreedt. 
d. De programmering moet op eenvoudige wijze in 

ALGOL verricht kunnen worden. 
e. Het geheel moet aangepast zijn aan de TR4 reken- 

machine in Delft. 
De in dit verslag beschreven simùlatiewijze, DIsAR 
genoemd, verschilt van de in de litteratuur ([1] t/m 
(6]) gevonden systemen voornamelijk in de interpretatie 
van wat hier onder punt b genoemd is. (Afgezien van de 
punten d en e, die verband houden met de in Delft ge- 
boden mogelijkheden) In genoemde litteratuur is het 
nauw aansluiten bij ae analoge rekenwijze alleen ver- 
wezenlijkt voor in- en uitvoer; meestal bestaat de 
invoer uit een lijst van elementen en hun verbindingen, 
de uitvoer uit grafieken of getallenreeksen van de 
gewenste uitgangen. ve berekeningen echter gaan (voor 
zover de schrijver dezes heeft kunnen nagaan) aan de 
hand van een stelsel uit de invoer gegenereerde dif- 
ferentiaalvergetijkingen. In de rekenfase van die 
programva's zijn de integratoren niet meer afzonder- 
lijk te onderscheiden. bovendien wordt in het hier 
beschreven systeem getracht de stuplengte te regelen 
aan de hand van de fouten in de integratoren afzonder- 


“lijk. Het enige systeem met staplengteregeling wat in 


de litteratuur gevonden is [9] heeft een staplengte- 
regeling die gebaseerd is op de totale fout van het 
systeem. 

Alvorens te komen tot de uiteindelijke vorm van net 
DISAR-systeem en het bespreken van de daarmee bereikte 
resultaten, zal eerst een eerder toegepast, eenvoudiger 
systeem besproken woraen omdat daarbij de moeilijk- 
heden op eenvoudiger wijze naar voren komen. 
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II. D& INTEGRATIEMETHODE VAN EULER. 


Van de elementen die een analoge rekenmachine kan 
omvatten zijn de integratoren de enige elementen 
waarbij de uitgangsgrootheid zowel afhankelijk is 
van de ingangsgrootheden en de tijd, als van de 
toestand* van het systeem (differentiatoren zijn 
buiten beschouwing gelaten). De overige elementen 
zijn elementen zonder toestand, de uitgangsgroot- 
heid daarvan is op een bepaald tijdstip bekend als 
de ingangsgrootheid (heden) op dat tijdstip gegeven 
is (zijn). Daarom is de integrator het enige element 
dat stappen neemt, Een stap van de integrator legt 
verband tussen de ingang y' op het tijdstip t en de 
uitgang y op het tijdstip t+h waarin h de staplengte 
is. 

Begonnen is met een zeer eenvoudige integratie 
methode, die aanneemt dat 

di y' dt=hy} (h=tj-t) zie fig. 1. 

Dit wordt de integratiemethode van kuler genoemd. 


Ô er 
Fig. 1. ĳ ee 
De integratie 4 opp.=h io 


methode van 
Euler. 
> é 


In de figuur is reeds te zien dat de fout evenredig 
zal zijn met y'' (zie III). Wanneer de formule 


Ey! 
Je Ae nd ho JI) gebruikt zou worden dan zou de 


en 
è 


gemaakte fout kleiner zijn (zie VI). Dit is ook in 


fis. l te zien. 


* Zie: Dr. Ir. H. Kwakernaak: 
Elementaire vegrippen uit de moderze systeemtheorie. 
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III. 


S= 


FOUTENBESCHOUWI 


Uitgaande van de formule van Taylor: 

ECtsh) = £(E) + hEI(E) + R (ea) 
met Re= h°(1-8) Ê"(t+H8) waarin O4GC1 

die gebruikt mag worden als f(t) twee maal gedifferen- 
tieerd kan worden, is de integratieformule 

F(t+h) = f(t) + hf'(t) (2) 
gebruikt. De fout die deze methode met zich mee 
brengt is dan gelijk aan R. Heeft de te integreren 
functie Zen “zlad verloop! dan mag aangenomen worden 
dat h°(1-6) r(t+eh) < nEe(t). 

In het vervolg wordt f'(t) (d.i. de ingang van de 
integrator) y(t) genoemd. Voor de fout per stap geldt 
dan: emhey'(t) . Over een interval a (bev. van tzt 
tot tzt ota) worden dan N =E stappen gedaan. De 
totale fout w is nu EN, hyt (t). Dus geldt: 


Eza.h.y! yr Ct) = a.h.y' (7) met LOST St ta (3) 


In het integratieproces wil men E niet te groot laten 
worden. Omdat echter in (5) de gemiddelde waarde van 
y'(t) niet bekend is zou de staplengte zo gekozen 
moeten worden dat "het redelijkerwijs altijd goed 
gaat'. Dit zou ecater onnodig veel rekentijd kosten. 
Daarom wordt in dit systeem een staplengteregeling 
gebruikt die zorgt dat E tussen een maximum en mini 
mum waarde blijft. Bit is bereikt door de staplengte 
h met % te vermenigvuldigen wanneer hy' Ct) groter is 
dan een bepaalde grenswaarde (max) en met 2 te ver- 
menigvuldigen wanneer hy'(t) kleiner is dan een 
andere grenswaarde (min). Bij elke stap moet er dus 
een maat aanwezig zijn voor hy'(t). Voor kleine h 
geldt: y(t+h) — y(t) = hy'(t), hetgeen direkt uit 
de definitie van de afgeleide volgt. Het kriterium 
(án de integratieprocedure is deze variabele 'krit' DE 
genoemd) is slechts te bepalen voor een reeds uit- « 
gevoerde stap, omdat y(t+h) de ingang van de inte- 
grator is op het tijdstip dat de stap t+h—>t+2h 
genomen kan worden. Dit is de reden dat de stap moet 
worden overgedaan wanneer blijkt dat hy'(t) max is. 
Hier is, evenals in deel IX terwille van de eenvoud 


niet de fout van het gehele systeem bekeken, slechts 
voor de integratoren afzonderlijk is een regeling 
getroffen. Het staat echter niet vast welke invloed 

dit op de fout zal hebben. De fout in het gehele systeem 
zal een functie zijn van de fouten in de afzonderlijke 
elementen en regeling van deze deelfouten zal de totale 


fout dus wel beperken. 
: 


IV. HET PROGRAMMA. 


Het ontworpen systeem, dat het simuleren van een 
door een analoog schema gegeven probleem mogelijk 
maakt werkt met de procedure's org en int. De pro- 
cedure int is de eigenlijke integrator, die als 
parameters o.a. de in- en uitgangsgrootheid (heden) 
heeft. De procedure int zorgt indien nodig tevens 
voor de staplengteverkleining (zie fig. 2b). De 
procedure org maakt voorbereidingen voor de volgende 
uit te voeren stap qn tevens voor de staplengtever- 
groting als dat nodig is. ken blokschema voor een 
gehele analoge opstelling geeft fig. 2a. 


Verbindingen van 
alle andere 
elementen onderling 
en met 
intesratoren. 


Beginvoorwaarden. 


Fig. 2a 


De opbouw van 
het digitale 
programma van 


integrator 1 


integrator 2 


een analoge 
schakeling. 


eventuele 
uitvoer 


De gebruiksaanwijzing van het systeem is als volg 


integrator n= 
integrator n 


l. Verklaring formele parameters van org en int. 
a. procedure int(in,uit,rc,k,p) 
Invoerparaneters: 
real in de ingangsgrootheid (integraná) 
real rc constante waarmee de ingangsgrootheid 
wordt vermenigvuldigd ( 1/7 in een 
analoge integrator). 
integer k rangnummer van de integrator 
label p maam van de label, die aan net begin 
van het "analoge!" programmagedeelte 
geplaatst moet worden (zie fig. 2b). 


Uitvoerparameter: 
real uit de uitgangsgrootheid. 
b. procedure org geen formele parameters. 
De opdracht org moet aan het einde van het 
“analoge! programma gegeven worden, vöór de 
uitvoeropdrachten. 
2, De beginwaarden. 
Hier volgt een lijst van identifiers, die door 
de integrator en de procedure org worden verklaard. . 
In de rest van het programma mogen dus geen vari 
abelen van deze naam worden gebruikt, 
store: een array waarvan de elementen als geheugen 
plaatsen worden gebruikt. 
t s de tijd. 
h : de staplengte. 
ho : de staplengte in de vorige stap. 
hmin : de minimale staplengte. 
min # minimum voor het kriterium. 
max ; maximum voor het kriterium. 
n s aantal integratoren. 


aa : 
bb : variabelen voor de organisatie van de 
dd H procedures. 


Vôór de label die aan het begin van het "analoge" 
programma staat moeten deze variabelen een waarde 
krijgen. De volgende variabelen zijn vrij te kiezen: 
hmin; de gewenste minimum staplengte. Het is aan te 
bevelen om h, ho en hmin dezeìfde waarde te 
geven. 
max en mins de te kiezen waarde van max hangt samen 
‘met de gewenste nauwkeurigheid. In de 
praktijk is gebleken dat de beste waarde 
van min ligt tussen %,max en 1/10,max. 
Beginvoorwaarden zijn in te voeren volgens: 
for iz=l step l until n do 
store [1,4] :=store [ì,9] t=bv, 
waarin bv, de i beginvoorwaarde is. 
De andere variabelen moeten een voorgeschreven waarde 
krijgen voor het begin van het "analoge! programma. 
De noodzakelijke opdrachten zijn: 


ij 


for is=l step l until n do 
store[i,1] :=store [Ì‚2] :=store [,3]:=0; 


ts=0; aas=bbsi=false; dd:=true; 


Fig. 5 geeft het analoge schema van een eerste orde 
systeem, bijlage l het gehele programma wat van 40'n 
systeem de responsie op een eenheidssprong simuleert, 
fig. 2b is het stroomdiagram van dit systeem. 
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Fig. j. Een eerste orde systeem. 
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RESULTATEN VOOR SYSTEMEN MET EEN INTEGRATOR. 


&. zerste orde systeem, resp 


ie op sprong. (tig. 5) 
Fig. # geeft de absolute fout voor verschillende 
waarden van het kriterium (staplengte variabel). 

Het kriterium blijkt hier een goede maat voor de 

fout te zijn, want verdubbeling van het kriterium 
geeft verdubbeling van de fout. Fig. 5 vergelijkt 

de fout die gemaakt wordt bij variabele staplengte 
met de fout die gemaakt wordt bij een konstante 
staplengte. Het gemiddeld aantal stappen per tijds- 
eenheid is in beide gevallen hetzelfde gekozen. 


È. Testprogramua voor staplensteregeling. 
Geintegreerd is de functie f(t)=6sin t +sin 3t 

In fig. 6a is het verloop van de functie en van de 
staplengte in dit geval weergegeven. Het is te zien 
dat de staplengte klein is als de helling van de 
functie groot is. In IIl is berekend dat de fout 
Ez=hy"(t). Staplengteregeling is dus op de juiste 
wijze geschied. In fig. 6b zien we de fout bij vari- 
abele staplengte, vergeleken met de fout bij konstante 
staplengte. (Het gemiddeld aantal stappen per tijds- 
eenheid is weer gelijk gekozen.) 


De absolute fout van het ee 
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verschillende waarden van krit. 
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Fig.l6a. Verlooy van de functie 6 sin t +!sin 3t Ix variabele staplengte 
ì 6b, 


staplengte 


Het foutenverloop van deze 
Functie 
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ME4R INTEGRATOREN, 


Aan de hand van de formule van Taylor (zie formule 
TE, CE) 


ËCt+h) = É(CU) + hÊ'(E) + R (4) 

is via 

£(t-h) = f(t) — hf!(t) + R' SS 
en door hierin t door t+h te vervangen, af, te leiden: 
Ê(t+h) = CE) + DE!CE+D) — R" 5) 
Uit (4) en (5) volgt: 

ECht+h) = f(t) + Bleie) + £'(te)} — Rt (6) 


Waarbij geldt: 

Rr = 1-0) fE“Ce+en) — e"Ct+(1-0)n)} 

In het algemeen zal (wanneer de functie een "glad" 
verloop heett) i'(t+en) — ["(t+(1-6)h) klein zijn, 
zodat 

R* 4 n(1-BIL"(t+Eh) 

deweze dat de fout R* kleiner is dan de in III af— 
geleide fout. De in het hiervoor beschreven systeem 
gemaakte fout kan dus gezien worden als een tijd- 
verschuiving over %h als f'(t) gebruikt wordt in 
plaats van HL CE) + £1(t+n)} « Op de ingang van de 
integrator is echter alleen f'(t) als grootheid aan- 
wezig. In geval van twee integratoren kan het effekt 
van deze verschuiving kleiner gemaakt worden door de 
tweede integrator als ingang f'(t+h) te geven. Bit is 
toegepast voor een sinusgenerator (fig.7). De ingang 
van de eerste integrator is f(t), als uitgang komt 
er Ee! (tan) .®) Wordt deze waarde als ingang voor de 
tweede integrator gebruikt (bij de stap van tijdstip 
t naar tijdstip t+h!) dan is in te zien dat de uitgang 
van deze integrator, dew.z. f(t+h), redelijk gecompen- 
seerd is voor genoemde tijdverschuiving. 


Ra seh £(E) 


fig. 7. -£(t) 

De sinusgenerator. 

%) Hier is verondersteld dat de integratoren met teken— 
omkering werken. Dit is bereikt door in alle formule's 
h door -h te vervangen. 


lj 


Voor het tweede orde systeem (fig. 8) is hetzelfde 
principe toe te passen, hoewel de compensatie niet 


zo goed zal zijn, omdat ook de uitgang van de eerste 


. integrator teruggevoerd wordt. 


fig. 8 
Het tweede orde systeem. 
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VII. RESULTATEN VOOR SYSTEMEN MET TWEE INTEGRATOREN, 


Voor een sinusgenerator is de uitkomst aan de uitgang 
van het systeem exact te berekenen. In bijlage 2A is 
dit gedaan voor de Euler methode, terwijl 2B de resul- 
taten geeft met de genoemde "tijds"correctie. De uit- 
komsten zijn: 

geval A: y' (uitgang 1 integrator) = 


(asn2)P/eE sin(t/h arctg h) (7) 


„& t/2h 5 

geval B: y' = eRe, 3) sine/n arctg(he} )} (8) 

In onderstaande tabel zijn amplitude-vergroting en 
faseverschil na 100 perioden uitgezet als functie van 
de staplengte voor beide gevallen, zoals gevonden 
uit de experimenten, Deze waarnemingen klovpen met de 


berekende uitkomsten. 


Afwijkingen na 100 perioden. 
Beginamplitude 1 frekwentie 1/2 per./sec. dus 


Be 7 5 stappen 
staplengte 1/10 betekent OH Gertode: 


fase- 
verschil 


7000000 
4,8 
1.17 


Methode A is minder nauwkeurig dan methode B. Een 
nevenverschijnsel wat methode B vertoont is echter 
dat de uitgang van de tweede integrator %h in fase 
voor loopt op de fase van de eerste integrator, ter- 
wijl bij methode A de fasen van eerste en tweede 
integrator gelijk zijn. Zowel in (7) als in (8) neemt 
de amplitude expomentieel met de tijd toe en het 
faseverschil lineair met de tijd. 


B. 


ede orde systeen. 


De berekening van de uitkomst is hier niet zo een- 
voudig te doen. De verkregen fout is in fig 9 aen b 


ide 


uitgezet voor verschillende staplengtes. In fig. 9b 
‚zien we evenals in fig. 4 dat verdubbeling van het 
kriterium een verdubbeling van de fout veroorzaakt. 
Voor dit systeem geldt evenals voor het in W.A be- 
sproken eerste orde systeem dat bij een konstante 


staplengte lin E = O waarin E de fouts Dit iseen 
t— co Bo chape 


eigenschap van elk systeem waarbij de uitgangswaarde 
naar een limietwaarde convergeert die onafhankelijk 
is van parameters als RC-tijd, demping en staplengte. 
Bij variabele staplengte zonder bovengrens zoals 
hier gebruikt, zal de staplengte zo geregeld worden 
dat het kriterium tussen twee grenzen wordt gehouden. 
Dat wil dan zeggen dat de fout nooit U wordt, omdat 
de afgeleide van het ingangssignaal boven een minimum 
wordt gehouden. 
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Fig. 9b. De absolute waarde van de absolute fout 


1 s van een tweede orde systeem. 
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VIII. DE PREDICTOR-CORRECTOR INTHGRATIEMETHODE, 


De tot dusver beschreven integratiemethode is niet 
alleen vrij onnauwkeurig, maar heeft bovendien het 
nadeel dat bij de gekorrigeerde methode een univer- 
sele programmering onmogelijk is, immers bij een 
willekeurig aantal integratoten moeten i integratoren 
als ingang jCt+h) en j integratoren als ingang y4(6) 
krijgen, Hoe deze integratoren verdeeld moeten worden 
is niet algemeen te bepalen. Deze nadelen kunnen 
worden opgeheven door een predictor-corrector methode 
te gebruiken, De nauwkeurigheid wordt dan bovendien 
groter, doordat in de correctorstap de invloeden 

van de verbindingen tussen de integratoren onderling 
merkbaar worden. Fig. 10 geeft een overzicht van de 
opbouw van een analoug programma gesimuleerd met een 
predictor=corrector methode. 


Verbindingen 
tussen alle 
elementen. 
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Fig. 10. Opbouw van een "analoog! programma 


voor een predictor - corrector methode. 


Het ge 
achte 


ozen integratiesysteem doet in een cyclus 


eenvolgens een predictor- en corrector-stap 
voor de halve staplengte, en met behulp daarvan een 


predictor- en correctorstap voor de gehele staplengte 


En 


Evenals het in de delen II t/m Vil beschreven systeem 
wordt de nauwkeurigheid geregeld door keuze van de 
stapgrootte. Per integratiecyclus worden dus steeds 
vier stappen gedaan, onafhankelijk van de nauwkeurig 
heid. 

Voor de afleiding van het invegratiesysteen wordt 
uitgegaan van de methode van Simpson: 


22 h EEn 
J adel XI? + + ELXos Ys zi 18e 
Je vraa} fl(zorvo) WEGEN) + £(2492)} (ate fie. 11) 


Hierin is f(x,y) = Y 


3 KX Has X2X0= 2h. 


de al 2 ax (tijdas) 
ed 
2h 


Fig. ll. Bij de integratiemethode van Simpson. 


‘In deze formule wordt verondersteld dat de waarden 
van de onafhankelijk variabele Kor Xj-CR Xp bekend 
0'lo)= Y5 3 fCjsrj)= YJ en 

fe, yod= Yh ; de waarden van de athankelijk vari= 


zijn, evenals f(x 


abele. Voor de integrator die hier bedoeld wordt 

mag de onafhankelijk variabele bekend worden veronder- 
steld,- immers HEX) X=XpXj waarin h de staplengte is. 
Verder zin aan het begin van de stap, dus in het 

punt XxX, de waaraen van y, en do bekend: 5 als uit 
gang en de als ingang van de integrator. 

De eerste opgave is nu het vinden van een schatting 
voor Ji Hiervoor doet de integrator de eerste stap 
(predictor): zit — YotLI5 ‚ waarin yi de voorspelde 
waarde is van Je Deze stap is identiek aan de stap 
van de integratiemethode vau kuler. Zoals in fig. 10 


oi 


is geschetst worden nu de verbindingen tussen alle 
elementen gelegd, waaruit nieuwe waarden voor de in 
gangen van de integratoren volgen. Deze ingangen worden 
hier Ji genoemd. Nu volgt de tweede stap (corrector): 
Y1Tot BOD: Dit is te vergelijken met de in 

deel VI onder (6) gevonden stap die in de literatuur 
[7] de trapeziumregel genoemd wordt, 

Uit de verbindingen tussen alle elementen volgen nu 
als nieuwe ingangen van de integratoren de waarden 

van Yi De waarde van Ya wordt nu gevonden door 
extrapolatie van een polynoom van de derde graad 

P‚(0), die voldoet aan: 

Pz) = Jo 5 Ps) = IJ t Pz) =Ij en PZG) =Jj't 
Voor deze polynoom geldt: 

Ps) = 5Jo#Itehyjt#hyj. In de tweede stap is 


gevonden: yj= Yot Bg) dus 
-2hyj'= HY+#I teh zodat de derde stap wordt: 


J = JoTeaIT'+#hyj . Uitv de verbindingen tussen de 
elementen wordt nu 5 gevonden waarna de vierde stap 
wordt gemaakt: Ya= zoenen ‚ overeenkomstig 

de methode van Simpson. 

Omdat de vier stappen hier beschreven "naar buiten! 
optreden als één stap, die in het vervolg cyclus 
genoemd wordt, ter onderscheiding van de stapgedeelte!'s, 
wordt in bovengenoemde formule's h vervangen door 
h/2, zodat de staplengte h de betekenis heeft van de 
afstand bussen Xx, eu Xoe Bovendien is de teken 
omkering en de tijdconsvante van een analoge inte- 
grator ingebracht door in het beschreven systeem 
voor h steeds -RC‚h in te vullen. RC heeft hier dan 
de betekenis van IM in een analoge integrator. 

De formule's zijn nu: 


see stap Yj = Vak 0.896 geeït als nieuwe ingangen vl 
2e stap Fr dT ROIS zi) taenn u u DE 
oid stap 5 So RC.n.(yj'-23j) nj Li ei ga 
4® stap Vas Tor ROB. Cypa tI5' laatste stap der cyclus. 


| 


De integratieformule's van het gebruikte systeem 
zijn terug te brengen op wat in [7] de eerste 
integratiemethode van Runge Kutta wordt genoemd. 
Noem Eje hy6 5 ko= hyj's Kz=Dy| en ky= hy! 
dan is de integratiemethode gegeven door: 

kj= hf(X09J0) 

kg= hf(xoth/2,Jotkj/2) 

kj+K 

kz= nie tb/2 te) 
- ky= hf +BrIoTkotêkz) 

E 1 
waarna volgt Jj= Jot ge täkztk,) 
Dit zijn de in [7] bij de eerste methode van Runge 
Kutta gegeven iutegratieformule's. 
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FOUTENBESCHOUWING EN STAPLENGTERKGELING. 


De nauwkeurigheid van het Runge Kutta integratie- 
proces is tot dusver slechts in bepaalde gevallen te 
schatten. Daarom is het eveneens voor het in deel VIII 
afgeleide systeem moeilijk een foutenbeschouwing te 
geven die aan de staplengteregeling ten grondslag 
ligt. In het volgende wordt echter getracht een voor 
de meeste gevallen redelijke benadering van de fout 

te vinden, die geleid heeft tot een stapleugteregeling 

die in de praktijk gunstige resultaten heeft opgeleverd. 

Allereerst volgt hier een lijst van veronderstellingen 

waarvan is uitgegaan. 

1e Er is verondersteld dat de fout die optreedt aan 
de uitgang van een gesimuleerd systeem beheerst 
wordt door van iedere integrator afzonderlijk 
de fout binnen bepaalde vastgestelde grenzen te 
houden. (Deze veronderstelling is reeds aan het 
slot van deel IIl genoemd.) 

2” &r is aangenomen dat de hogere afgeleiden van de 
functie die beschouwd wordt convergeren, en wel zò 
dat HEEN CPE) woor arte nne overs 
deel VI gestelde eis dat de functie een "glad" 
verloop heeft komt in de praktijk op hetzelfde neer. 
De voorwaarden moeten gelden vocr elke t tussen 
het begin en eindtijdstip van de periode waarover 
geintegreerd wordt. 


NN 


Er is verondersteld aat de fout in elke integrator 
ongeveer gelijk is aan de afkepfout van de Simpson 
integratiemethode hoewel de gebruikte methode 
slechts geschatte waarden van Yi en y5 invult in 
de formule van Simpson. 

De afkapfout bij de simpsonmethode is 


MCS) E hs zh 8 EE > 
B = 50 ETT) met 14 T4t+h (als veronderstelling 2 
geldt!). De in een tijdsinterval T van to tot t+T 
gemaakte cumulatieve fout ig dan 


ie 
En 50 £ A) met CT, tussen ten tj Zodat 


BE 
Ez 50 LATE) met € tussen toent rr (9) 
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Getracht wordt de fout onder controle te houden door 
de staplengte te regelen. In elke Sels van vier 


5 stappen wordt Seed eh of maxph ae >min waarin 
5 AÊrantn de tweede laifferentie van het ingangssignaal, 
—_ n en max en min van te voren vastgestelde grenzen zijn; 
z die een maat voor de nauwkeurigheid zijn. Omdat de 
E, fout evenre is met Pr) zou men kunnen verwachten 
Rd dat voor deze test een hogere differentie gebruikt zou 
moeten worden. Er zijn vijf redenen te geven waarom 
p dit niet is gedaan. 
] 1. Af is zonder extra rekenwerk bij de vierde stap 
ai der cyclus te verkrijgen, n.l. SE zi VoT2IFTË' « 
@. Het berekenen van hogere differenties gaat met 
) extra fouten gepaard, waardoor deze hogere dif- 
=S ferenties niet altijd een betere maat voor de fout 
E zullen zijn. 
5. De veronderstelling dat de fout in dezelfde orde 
] van grootte zal liggen als de fout bij een Simpson 
= integratie methode ís alleen dan waar als de schat 
ting van Yi goed is (eerste en tweede stap). 
ad fout in deze schatting is ongevecr evenredig met 


f'(T). Dit is in te zien door te vergelijken met 
de in deel VI gevonden formule (6). 
4e Het is aannemelijk dat de nauwkeurigheid mede in 
z e Eee 
verband zal staan met de kromming WO d 
5e Onder de veronderstelling dat rede) 
geldt dat bij regeling aan de hand van de tweede 


jm je 


differentie de staplengte waarschijnlijk kleiner 


5 


dan nodig wordt gehouden en dat de nauwkeurigheid 


l 


dus niet kleiner zal zijn dan bij gebruik van de 

_ vierde differentie. 

In deel VIII is nog genoemd dat voor nh overal -RC.h 

is ingevuld, om een goede overeenkomst te krijgen met 

analoge integratoren. Eet foutenkriterium wordt nu 
RO) NELE. Dit foutenkriteriun, met een verdubbeling 


van de staplengte als (RC) PhEASE Ê4min voor elke 
eeen en halvering van de staplenste zodra in 
ee ên der integratoren (£ C)nSAs f>max geeft bevredigende 


AB 


resultaten, Hierbij was max een door de gebruiker te 
kiezen invoergrootheid, min werd steeds max/40 genomen. 
In het gebruik is gebleken dat de maximale fout in de 
grootte orde van de gekozen waarde van max blijft en 
dat de maximale fout van een bepaald systeem bij 
verandering van max evenredig is met de gekozen waarde 
van max. Vergelijk de resultaten beschreven in deel XI. 


X. EIGENSCHAPEEN VAN HET DISAR VOORPROGRAMMA. 

De hier ontwikkelde integratieprocedure is verkrijg- 

baar in een voor programmering in ALGOL op de TR4 

algemeen bruikbare vorm; de vorm van een voorprogram- 
ma, leverbaar op een los telex-code ponsbandje. 

Allereerst volgt hier een korte verklaring van de 

ALGOL-structuur. Het voorprogramma (bijlage 5) bestaat 

uit een niet afgesloten "blok" en een niet afgesloten 

binnenblok''. In deze twee blokken zijn de "verklaringen" 

(declarations) van de variabelen (identifiers) en de 

procedure opgenomen. Daarna volgt een aantal opdrachten 

die de verklaarde variabelen een waarde geven. Gevolgen 
van deze structuur zijn: 

a) Alle door de gebruiker verklaarde variabelen zijn 
lokaal ten opzichte van de in de integratoren ge- 
bruikte variabelen. (zie eigenschap 1la) 

b) In het programma van de gebruiker moeten deze twee 
blokken afgesloten worden, dus is twee maal extra 
het woord end verlangd. 

ce) Voor het programma van de gebruiker mag geen stuur- 
informatie (8) staan; de TR4 moet het programma 
van de gebruiker opvatten als het vervolg op het 
DISAR voorprogramua. 

Het DISAR voorprogramzua verlangt n+5 getallen op de 

‚getalband (vóêr de eventuele door het programna van 

de gebruiker verlangde getallen). Hierin is n het 

aantal integratoren. De verlangde getallen zijn: 

LS getal is het aantal integratoren dat wordt gebruikt. 

2° getal is de gewenste minimum staplengte (zie de 

eigenschappen 2, 5 en 4). 
5° getal is de gewenste waarde van de variabele max 
(max as de maximale absolute fout). 

uS t/m Dit zijn de beginvoorwaarden voor de integratoren 

n+3° getal 1 t/mn d.w.z. de uitgangswaarden op 

het tijdstip t=0. 

In fig 10 is reeds de opbouw van een programma van een 

gesimuleerd analoog schema geschetst. De getekende 

integrator heeft 2 stapgen in de integratiecyclus, 

de in DISAR gebruikte integrator heeft een cyclus 


van vier stappen (zie fig. 12). In fig. 15 is geschetst 
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Ja 
bereken 


Ee bewaar JT bewaar Ji nauwkeurigheid 
T krits=yl-eyl +5 
Q 


foet de stap 
i= Jor verkleind ? 
5 ee, 5 5 ot a nde 
J18=Jo 7 hyj'+2hyj 


geef aan als 
krit te klein 


verklein 
nns= 15 


alle uitge 
1ieuwe wa: 


alle uitgangen 
nieuwe waarden) 


vergroot 


uitgang 


Fig. lë. Blokschema van een integrator. 
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Declarations en statements 
teb.v. het programma van 
de gebruiker. 


il “vast'' programma 
gedeelte. 


Verbindingen van uitgangen 
der integratoren, via overige LE 
elementen naar de ingangen 


“repeterend'' 
Integratie opdrachten. gar 


programma gedeelte. 


Uitvoer opdrachten 


Sprongopdracht naar label. Vv 


nse cijfers korresponderen met de 
cijfers in bijlage 4 t/m 6. 


hoe het programna van de gebruiker er zal komen uit 
te zien. áls voorbeelden mogen verder dienen de bij- 
lagen 4 t/m 6. 
EIGENSCHAP 1. zen integrator wordt gesimuleerd door 
de procedure opdracht int(in,rc,k,lab) 
a) Verklaring formele parameters: 
real in Dit is de invoergrootheid. Als actuele 
parameter moet een getal (number) of uitdrukking 
(expression) worden ingevuld, die de ingangs- 
grootheid van de integrator simuleert. 


real re Invullen de konstante urmee bovenstaande 
parameter vermenigvuläigd wordt. (rc simuleert 
1/C van een analoge integrator.) 

integer k Invullen het rangnummer van de integrator. 


label lab Invullen de door de gebruiker gekozen 


iabel die Oeds lí is aange=- 


naam van de 
geven. 
b) Uitvoerparaneter. 


Het array element uit[k] geeft de getalwaarde 


2305 


van de uitgang van de ke integrator aan. Deze is 
te gebruiken in deien II en IV van fig. 15 en 
bovendien in de 'exvression' die voor de onder a) 
genoemde invoergrootheid wordt ingevuld. 
Uit bijlage 5 is af te leiden dat de bepaling van de 
tijd t geschiedt door na elke cyclus t met h te ver- 
hogen met de ALGOL opdracht ts=tt+h; Omdat getallen in 
het algemeen bij de interne representatie in de TR4 
als een (binaire) repeterende breuk verschijnen en 
het aantal bits van een woord in de T&4# constant is 
ligt de fout bij elke optelling tussen 0,013.1071% 
en 5.10"2% 
zo'n geval kan bij integreren van t=0 tot t=256 met 
een staplengte van O.Cl de absolute fout in t tot dent 


van de waarde van het grootste getal. In 


toenemen. Dit is geconstateerd bij de meting van de 
faseverschuiving van een sinusgenerator. De fasever- 
schuiving was in sommige gevallen kleiner dan de fout 
in de bepaling van t. Wanneer echter een getal voor 
de staplengte gekozen wordt dat in de interne repre- 
sentatie van de TR4 een echte (niet repeterende) 
binaire breuk is zal de fout tot O gereduceerd kunnen 
worden. 

Wanneer nu de gebruiker niet zelf door het invoeren 
van een minimum staplengte de beginwaarde van de 


staplengte bepaalt (zie eigenschap 2), dan wordt de 


s ae eerste 5 signi= 


ficante binaire bits van nul verschillen. Dit gebeurt 


met de opdrachten: 


entier(tzefnn) while hg & do nn:=nn+l; 


hs=he2f( on); 
hmins 


(hier zijn t en nn als hulpvariabelen gebruikt.) 


Nu zijn de drie volgende eigenschappen van DISAR bekend. 


EIGENSCHAP @. wanneer de gebruiker als minimum stap- 
lengte op de getalband een waarde kiest die groter 


ie) 


is dan 10 ° dan wordt deze waarde als beginwaarde 


van de staplengte gekozen. 
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EIGENSCHAP 5. Wanneer de gebruiker (kleine) afrondings- 
fouten in de bepaling van de tijd t wil vermijden, 
dan zal hij als minimum staplengte moeten invoeren 
of een waarde kleiner dan 10° (bve O0) of een 
echte binaire breuk b.v. 2° waarin n geheel en 
(meestal) negatief, 

EIGENSCHAP 4, Wanneer de gebruiker als minimum stap 
lengte 0 invoert, dan wordt een niet-repeterende 
breuk als beginwaarde voor de staplengte gekozen. 
Door een discontinuiteit of iets dergelijks in een 
der ingangsfuncties van de integratoren, of door 
een te klein gekozen waarde van max kan het in dit 
geval voorkomen dat de staplengte steeds verder 

.verkleinà wordt zonder ooit een bevredigende waarde 
te vinden, Voor de meeste systemen echter is dit 
gevaar niet aanwez 


Het mechanisme wat zorgt voor het vergroten van de 
staplengte bestaat uit de opdracht aas=(krit{minjAaa 
âie in elke integrator wordt uitgevoerd, en de opdracht 
if aa then hs=2sh die alleen door de laatste integrator 
wordt uitgevoerd. Hieruit volgen eigenschap 5 en 6. 
EIGENSCHAF 5. Wil de gebruiker dat de staplengte 

gooit vergroot wordt dan kan hij in het repeterende 

deel van het programma de opdracht aa:s= false geven. 
EIGENSCHAP 6, Wil de gebruiker de staplengte aan een 

meximum waarde binden dan kan hij als eerste opdracht 

na de label aan het begin van het repeterende 

programmagedeelte de opdracht plaatsen: 

ZE hjx then 


aa:= false end; 
Hierin is x de gewenste maximum staplengte. 
Afgezien van de in te vullen getallen op de getalband 
behoeft de gebruiker geen voorzieningen meer te treffen 
die de organisatie van ce integratieprocedure behelzen. 
Het is echter wel mogelijk alle variabelen naar eigen 
keuze in te voeren. 
EIGENSCHAP 7, Wil de gebruiker een beginwaarde van een 
variabele zelf kiezen, dan kan dit gebeuren in het 


vaste programmagedeelte (voor de label). 


a). De waarde van t aan het begin kan worden gekozen 
door de opdracht ts=tè:za waarin a een willekeurig 
getal (eventueel negatief). 

b). wen beginvoorwaarde voor een der integratoren 
kan worden veranderd door de opdracht 
store [p‚2]:=uitfp]:=bv; waarin p het nummer van 
de betreffende integrator en bv de gewenste begin- 
voorwaarde. 
ce). Alle andere variabelen kunnen onafhankelijk 
van elkaar worden gekozen. 


Enkele eigenschappen die verder geen verklaring be- 

hoeven zijn: 

EIGENSCHAP &, De staplengte kan ook ester worden 
gekozen, echter alleen bij een constante staplengte. 
Wordt in het vaste programmagedeelte de opdracht 
hs=hmin:=a gegeven met a negatief, dan is a de 
gewenste constante staplengte. Gevolg is dat de 

ze t wordt steeds kleiner. 

EIGENSCHAP 9. Loor invullen van een negatief getal voor 


“tijd terugloopt'', de 


RG in de procedure-ovdracht int wordt de teken 
omkering opgeheven. 

EIGENSCHAP 10. De waarde van t kan (ook) in het repe- 
terende programuagedeelte door de gebruiker wille- 


keurig veranderd worden zonder invloed op de 
werking van de integratieprocedure als zodanig. 
EIGENSCHAP 11. Afgezien van een variabele met de naam 
“uit! kan de gebruiker variabelen verklaren met een 
willekeurige naam. Dus ook namen als h, t, max 
en dergelijke kunnen worden verklaard en gebruikt 
zonder invloed op de werking van de integrator. 
EIGENSCHAP 12, De simulatietijd van een systeem 
(rekentijd van het programm ) hangt af van de 


staplengte, het aantal integratoren, de rekentijd 
voor de overige ejementen e.d. slechts onder 
voorbehoud kan nier medegedeelc worden dat een 
integratiecyclus (van vier stappen) per integrator 
0,210 ee min. duurt, Voor een systeem van vijf 
integratoren is dus waarschijnlijk een snelheid 
van 1000 stappen per minuut te halen. 
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RWSULTaTeN BEREIKT MeT HeT DISAR VOOREFROGRANMA. 


Niet alleen ter vergelijking met het in VII A genoemde 
resultaat, maar ook als op zichzelf staande kwaliteits 
beoordeling is de sinusgenerator geschikt. 

Onderstaande tabel geeft een beeld van de amplitude- 

en faseafwijkingen na 100 perioden bij verschillende 
staplengtes (beginamplitude 1, frekwentie 1/2M peri- 
oden per seconde). Hierbij moet rekening gehouden 
worden met het feit dat niet de frekwentie, maar het 
aantal stapven per periode de nauwkeurigheid bepaalt, 


en dat de staplengteregelin; 
aantal sta 


t dat het 


r voor zorg 


ppen per periode bij verschillende frekwenties 
bij dezelfde waarde van max steeds minder dan een factor 


„twee onderling atwijken. 


staplengte taseverschil. 


in radialen. in perioden. 


niet meetbaar* 
=5.0 106 [eo 207 
niet meetbaar* 
=5.0 10 =8.0/ 1072 
elen a nhene 


Het verloop van de afwijkingen met-de tijd vinden w 

in fig. 14 a en b. Vergelijken we deze grafieken met 
de in VII A gevonden resultaten (6) en (7) dan valt 
niet alleen de grotere nauwkeurigheid op, maar boven- 
dien een fout die niet meer exponentieel, maar lineair 
toeneemt met de tijd. Dit is in overeenkomst met het 
in IX formule (9) berekende foutenverloop. 

Bijlage 4 geeft het aLGOLprogramma. 

: 3e variabele gex op de nauwkeurigheid. 
In fig. lva is de fout in amplitude en fase, bij de 

* Wegens ongeschikte keuze van de steplengte, 
zie deel X, eigenschap 3. 
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sinusgenerator na 100 perioden, uitgezet als functie 
van de ingevoerde waarde van max. in dezelfde grafiek 
is het verloop van de maximale fout van de responsie 
van een gesimuleerd tweede orde systeem op een stap- 
functie uitgezet.(zie ook XI D) De waarde van max blijkt 
een redelijke maat te zijn voor de grootte-orde waarin 
de fout zal vallen. Het verband tussen de fout en de 
ingevoerde waarde van max is echter niet eenduidige 
(vergelijk ook met XI F) Voor andere systemen zijn 

de fouten dan ook nog niet in hun afhankelijkheid van 
max te voorspellen, maar moeten zo mogelijk experi 
menteel worden bepaald. Omdat in het algemeen de 
uitgangsspanning van een analoog systeem niet is te 
berekenen, is het maken van een dergelijke grafiek 
uiterst moeilijk voor een willekeurig systeem. 

Fig 15b laat een ander effect zien. Hier is voor een 
sinusgenerator de afwijking na 10C perioden als functie 
van de gekozen beginwaarde van de staplengte uitgezet, 
bij een constante waarde van max. De staplengteregeling 
houdt de nauwkeurigheid bij max = 1073 tussen 10 en 
4.102. Omdat de staplengte discontinu geregeld wordt 
en wel door vermenigvuldigen met 2 of met % is deze 
variatie vrij groot. Om het bereik van de staplengte 


‘ groot te houden is een fijnere regeling niet wenselijk. 


Y. Iavloed van rc en van de grootte van het ingangssignaal. 


De in fig. 16 getekende analoge schakelingen geven 
dezelfde uitgangsspanning als Ì 


schakeling a de 
ingangsspanning Asin wt is en bij b lO.Asin wt. In de 
integratoren van het DISAR-systeem ís ook een tijd- 
constante in te voeren en hetzelfde effect trecdt daar 


op (wanneer rc gemaakt wordt aan 1/Y = 1/RC 


en Et 
OISAR zet 
resp. 1O/RC). Bij net DISAR-systeem treedt echter de 


Verschillende ingangsspanningen die dezelfde 
uitgangsspanning veroorzaken. 
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Fig. 15a. Fouten als functie van 


Systemen. 
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bijzonderheid op, dat de nauwkeurigheid in beide 
gevallen niet aezelfde is. Het kriterium voor de 
stavlengteregeling is namelijk (ro)Pn2a' Se waarin 
At de Ges differentie van het ingangssignaal; 
dus (re) on ENE zee) onten, Wordt nu PET SAR N maal 


zo groot dan wordt ook £' N-maal zo groot, zodat 


het kriteriwa No maal zo groot wordt. Wordt echter 
het ingangssignaal met N vermenigvuldigd, dan:heeft 
dit slechts effect op f'' zodat het kriterium slechts 
N maal zo groot wordt. De staplengte wordt in beide 
gevallen niet op dezelfde wijze geregeld, hetgeen, 
zoals in fig. 17 te zien is, leidt tot een verschil 
in nauwkeurigheid. Het is dan ook uit oogpunt van 

de nauwkeurigheid aan te bevelen om de waarde van 
re zo te kiezen dat de ingangssignalen niet extreem 
groot worden. Van "oversturing" zoals in het analoge 
geval is hier echter geen sprake. 


D. Responsie van een tweede orde systeem. 


In fig. 18 zijn absolute en relatieve fout uitg 
voor de responsie van een kritisch gedempt tweede 
orde systeem op de eenheidssprong. De waarde van 
het kebblium was 2,146.10”/. Opvallend bij dit systeem 


was dat de vorm van de foutencurve voor de verschillende 


waarden van max bijna identiek was. Dit is als volgt 
te verklaren: De gemaakte fout hangt in hoge mate af 
van de staplengte, de vorm van de toutenkromme hangt 
dus af van de punten waarop de staplengte verandert, 
Nu is voor elke volgende meting ae beginwaarde van h 
half zo groot gemaakt en max 1/16 maal zo groot. 
Wordt nu nieu vergeleken met max, dan overschrijdt 


het kriterium steeds in hetze punt de waarde 
vau max. De gevonden foutencurven moeten dus inderdaad 


vergelijkbaar zijn en bovendien is nu uit het verloop 


van de verhouding maximale fout / waarde van max 


te zien of hier inderdaad sprake is van een nauw= 
EN N zi > 5 

keurigheid evenredig met h'. In onderstaande tabel 

waarvan de eerste twee kolommen reeds voor een lijn 


in fig. 15a zijn gebruikt, is te zien dat de ver— 
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Fig. 18. Fouten bij responsie van een tweede orde systeem op sprong. (max = Zele: 16-73 
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Verder is een tweede orde systeem met een kwaliteits- 
factor Q=35.5 en een natuurlijke frekwentie wj=225 Hz 
(zie fig. 19) gesimuleerd. De frekwentiekarakterestiek 
van het systeem is bepaald, waarbij de gevonden 
maximum amplitude en bijbehorende frekwentie over- 

een kwamen met de berekende waarden. Evenals bij een 
dergelijk analoog systeem is gevonden dat na iedere 


frekwentieverandering een inslingerverschijnsel op- 
treedt, dat pas na een aantal perioden (plm. 5) is 
uitgedempt. Voor een echte analoge opstelling is deze 
invloed minder dan digitaal is-gevonden omdat zowel 
de analoge sinusgenerator als de analoge schakeling 


zoals in fig 19 getekend dempingen door parasitaire 


capaciteiten, lekweerstanden e.d. bezitten, die 


zeer hoge frekwenties uitfilteren. Een digitale 


frekwentieverandering veroorzaakt een duidelijke 


k in het ingangssignaal. Omdat steeds op het uit 
dempen van het inslingerverschijnsel gewacht moest 


worâen, duurde een bepa 


van twintig punten van 
de frekwentiekarakterestiek bij 1% nauwkeurigheid 
vier minuten rekentijd. Daarom zijn nauwkeuriger 


metingen achterwege gelaten. 


Het analoge scaema van de door dhr. G. Broekstra 


ontwikkelde opstel is geschetst in fig. 20. 
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Fig. 19. Tweede orde systeem. Responsie voor ver- 


schillende frekwenties. 


derde orde systeem — 


Beslis 
aha ser 
rid 


Fig. 20, Minimum zoekend systeem (een dimensionaal). 


Een beslisser moet de ingang van een functiegenerator 


zo regelen dat de uitgang van het derde orde systeem 


minimaal wordt, na een willekeurige beginvoorwaarde. 


De beslisser beslist aan de hand van de afgeleide 
van de uitgang wanneer de schakelaar bediend wordt. 
Een wachttijd is ingesteld zo dat âe schakelaar 
slechts met bepaalde tussenpozen kan worden omgezet. 
n de hand van de stand van de schakelaar wordt 


Aa 
+16 V of -16 V 


ntegreerd in de eerste integrator. 


Als functiegeneratoren zijn een absolute waarde nemer 


en een functiegenerator die de ingangsspanning kwa 


drateert gebruikt. De resulterende spanning diende 


s ingang van het derde orde systeem. Toen aan de 
ingang van de beslisser ruis werd toegevoerd bleek 
dat door het reeds in XI D (blz 43) genoemde effect 
de uitwerking van de ruis ernstiger was dan in het 
analoge 


val, mede omdat bij deze simulatie de 


beslisser traagheidsloos werkt. Het systeem heeft 
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=45- 


ook gewerkt in het drie dimensionale geval, d.w.z. 
het geval dat de beslisser drie ingangen KjsXg en Xs 
bedient, als waarde voor max is ingevoerd 107, 

als minimum staplengte is 10”“gekozen. Kleinere 
waarden van max en van de minimum slaplengte gaven 
verschillen die kleiner waren dan de door dhr. Broek- 
stra gewenste nauwkeurigheid. Het is tot op heden 
niet gelukt een exact verloop van de uitgang y theo- 
retisch af te leiden; de uitkomst is vergeleken met 
de uitkomsten van de analoge opstelling. Bijlage 5 
geeft het AIGOLprogramma van het drie dimensionale 
systeem. Meer resultaten staan vermeld in het kan- 
didaatsverslag van dhr. G. Broekstra, d.d. 6-4-'65, 


we 


E05. 


Een door Ir, H.G, Stassen gebruikte analoge opstelling 


“voor het simuleren van de bewegingen van het onder- 


stel van een treinwagon is gesimuleerá. Het schema 

bevatte: 

4 integratoren, 

8 fFunctiegeneratoren, waarvan 4 door een vergelijking 
gegeven (absolute waarde nemers) en 4 door een 
grafiek gegenen, 

l tijdsvertraging, 

4 deelschakelingen voor het delen van twee analoge 
signalen, 

2 relais, 

15 optellers en tekenomkeerders, 

25 potentiometers. 

De responsie op een stapfunctie is in fig. 21l ge- 

tekend. Logarithmisch decrement, nuldoorgangen en 

toppen waren dezelfde als gevonden met behulp van 

de analoge rekenmachine. Deze vergelijking levert 

echter geen informatie over nauwkeurigheden beter 

dan 1%. Daarom is getracht een indruk van de gemaakte 

fout te geven door resultaten behaald met verschillende 

waarden van max met elkaar te vergelijken. Zo werd 

de volgende tabel samengesteld. De in deze tabel 

uitgezette punten zijn bepaald door derde orde inter- 


enb 


tijdstip |tijdstip |tijästip |tijdstip 
ze 


1e nulpunt De nulpunt | 5” nulpunt 4enulpunt 


107%} 0.2007133 | 0.57005041 
107*| 0.2007231 | 0,57005623 
2| 02007235 | 0.57005706 

0.2007237 | 0.57005718 


1.508722 
1.508724 


tijdstip uitslag 


1° minimum 1° maximum 


0, 55669687 | 0, 91850054 | C,72601296 
0. 35670064 | C,91825153 | 0.726035546 [C,4199927 

„55670179 | C,91825989 | 0.72605535 |O,4199907 
0,55760187 |0.,91825912 | 0,72605535 [C,4199902 


Q.4200050 


polatie tussen de vier dichtstbijliggende berekende 


punten. Interpolatie tussen de punten a, b, c en d 
in fig. 22 van de 
afgeleide van de 
uitgang 
tijdstip 
het maximum optreedt, 


zeeft het 


vaarop 
interpolatie tussen 
de punten a',‚ b', c' 
en d' geeft de bij- 
behorende waarde van 
de uitgang. (Voor de 


bepaling van het tijd- 
stip van een nul- 


Fig. 22, Zie tekst. 


doorgang is slechts 

een interpolatie noâig). De staplengteregeling werd 
in dit geval geheel door DISAR verzorgd; als minimum 
staplengte is O ingevoerd. In fig. 25 is het verloop 
van de staplengte getekend voor de verschillende 
waarden van max. Het bleek dat de staplengte bij 

max = 10 en 1077 gaat hunten. Dit wordt veroorzaakt 


door de door het DISAR programma gekozen verhouding 
max 
Zn 
min 
doen verdwijnen. Een essencieel bezwaar ís dit hunten 


= 40, Een groter gekozen waarde zou dit effect 


echter niet. Bijlage 6 geeft het ALGOLprogramma. 
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tie treinbeweging. 


gteverloop bij simula 


25. Stapleng 


Fig. 
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G. Enkele anders systenen. 


Er is gesimuleerd hoe het gedrag van een eerste 

orde systeem zal zijn op het in fig. 24 getekende 
ingangssignaal. Afhankelijk van de tijdconstante 

van het eerste orde systeem loopt het (berekende) 
uitgangssignaal achter op het ingangssignaal en dempt 
het verschijnsel minder snel uit. Voor een ingangs- 


signaal zeen 


voor Tt, en O voor Dt, heeft de res- 
ponsie van een eerste orde systeem met £=1/20 een 
symetrische vorm. Een aardig bijverschijnsel is dat 

de tweede afgeleide van de ingang van de integrator 
continu is en dat dus de staplengteregeling het 
breekpunt "niet ziet'. Fig. 25 tracht een indruk van 
de nauwkeurigheid te geven. Het snijpunt van de twee 
rechte lijnen moet een meetvunt zijn, als het systeem 
goed gewerkt heeft. Dit was tot op hoge benadering het 
geval. 


— Ingangssignaal 


== responsie 


Fig. 24. Een gebruikt ingangssignaal voor een 
eerste orde systeem, met responsie. 


Tot slot kan worden medegedeeld dat Dr. Ir. H. Kwaker- 
naak het DISAR programma heeft gebruikt voor de op- 
lossing van de vergelijking van Ricatti: 


Ee EEn 


Hierin is Q de gevraagde 5x5 matrix, d een gegeven 


B'} + a + A'Q 


vector (3 componenten) en A, B en F gegeven matrices. 
Van belang waren de stationaire oplossingen. 
Het resultaat was bevredigend. 


iesponsie 
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XII, CONCLUSIES. 


Hoewel vooral deel IX veronderstelling 2 (blz. 26) 
een beperking betekent voor de klasse van functies 
die geschikt zijn voor het gebruik van de hier ge- 
schetste integratiemethode, is in de praktijk ge- 
bleken dat DISAR vrij algemeen bruikbaar is. 

Een andere beperking voor het gebruik is de rekentijd 
die sommige programma's eisen. Als vuistregel kan 
gelden dat het meestal niet zinvol is om een periodiek 
of bijna periodiek verschijnsel, waarvan meer dan 
honderd perioden onderzocht moeten worden, te simu- 
leren. De rekentijd hangt ook af van de gewenste 
nauwkeurigheid. Immers bij een grotere staplengte 

is weliswaar de rekentijd korter, maar de nauwkeurig- 
heid kleiner. Een (extreem) grote nauwkeurigheid 

is meestal niet erg zinvol, maar in sommige gevallen 
opent ze de mogelijkheid tot digitale verwerking 

van de analoge resultaten. 

Onder de uitbreidingen die mogelijk zijn kunnen 
genoemd worden het simuleren van meer elementen dan 
alleen integratoren, waardoor de voor de gebruiker 
benodigde kennis van ALGOL kleiner zou kunnen zijn. 
Het is echter mogelijk dat de flexibiliteit dan 
afneemt. De in deel XI E besproken opstelling, die 
gebruik maakt van logische bewerkingen en besturing, 
demonstreert deze flexibiliteit. 

Een andere uitbreiding zou zijns het automatisch 
schaalfactoriseren van een gegeven analoog schema. 
Bij de eigenschappen van de staplengteregeling zijn 
voor- en nadelen nog sterk verweven. Tegenover het 
voordeel van automatische staplengteregeling aan de 
hand van de grootte van max staat het nadeel dat het 
verband tussen de fout van de simulatie en max niet 
eenduidig ise 

Misschien mogen hier tot slot de woorden herhaald 
worden, die Dr. J.A. Zaat uitsprak op de openbare 


‘les gegeven bij de aanvaarding van het ambt van 


Lector in de nummerieke wiskunde aan de Technische 
Hogeschool te Delft, op vrijdag 1l december 1964: 


Vd ed md ml md el hd hd ed eed eed beed Led Ld Ld Li Lt La Lj 


EOD 


“Samenvattend kan men zeggen, dat op het gebied 
van de constructieve numerieke analyse en van de 
numerieke oplossingsmethoden in het algemeen nog 
veel werk gedaan kan worden. Ik denk in het bijzonder 
aan het convergentie-onderzoek dat met behulp van 
rekenmachines zou kunnen worden uitgevoerd, aan 
rekenprogramma's die de machine laten werken met 
zelf-corrigerende stapgrootte, om daarmee de sta- 
biliteit van het rekenproces te verzekeren, aan het 
grote belang van de beheersing van de foutenvoort- 
planting." 


mn nd mend 


] 
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zel and Jel. Wilson: 
Digital computor programs, analog solutions. 
Control Eng. June '6l page 155 


and J.J. Skiles: 


itally simulatea analog computor. 
ference 1965 


Froceedi Spring Joint Computor Con 


RA. Gaskill, Jeu. Harris and Al. Mc Koight: 


DAS: A digital analog simulator. 


Proceedings Spring Joint Computor Conference 25 1965 


1d bet. Warshawsky: 
al computation. 


An enlarged version of MIDAS. Simulation Okt. 1964 
R,D. Brennan: 


FPactolus:s A mulation language. 


jimulat 


Prof. D 


Collegedictaat a 15; Numerieke analyse. 


Fonsvoorschriften en uitvoer van ALGOLprogranma'!s 


(Stencil iskunaige Dienst d.d. 8 juli '64) 
H.E. Petersen, F.J. Sansom, e.a. 
MIDAS. Proceedings AFIFS 26 deel 1 (1964), 515. 
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« Programma met procedures dat de responsie 


op een stap 


n een eerste orde sys 
functie simuleert. 


Ì 
IBEGIN? REAL? T‚H,HO,HMIN,NIN,MAX; 
)ARRAY' STORE[1:2,1:61; 
INTEGER” N; 
*BOOLEAN” AÀ, BB, DD; 
— PROCEDURE” ÓRG; 
1BEGIN® "INTEGER’ 1; 


STORE [K,‚3]:=STORE[K, 41 =STORELK, VIER Gej 
UIT: =STORELK, 31 


8 'FOR' Ir=1 ÍSTEP! 1 UNTIL” N ’DO’ 
3 YBEGIN! STORELI,212=STORELI, 13; 
8 STORELI , 512=STORELI, 41; 
8m STORELI ,412=STORELI , 32 
5 PEND’; 
e Ts=T+H3 HO z=H; 
zo EE f 4 
E 1BEGIN? IE AA TEN 
E IBEGIN! Hr=2=H; BB2="TRUE';  AAs= 'FALSE'; 
5 PRINT(**H WORDT **,H) 
Ei 1END” JELSE’ AA z=/TRUE? 
2 1END’ 1ELSE' DD:= 'TRUE' 
1END’ 
B T1PROCEDURE! INT(IN,UIT,RC,K‚P); 
IREAL? IN,UIT,RC; 
"INTEGER? K; 
LABEL! P; 
1BEGIN? 'REÁL' KRIT; 
EN DD ZEN? 
STORELK, 1 z=1N; 
KRITz=ABS(STORELK, 11-STORELK, 21); 
'IE* BB *THEN’ 
TBEGIN! *IF? KRIT ”GREATER' MAX "THEN 
el 1BEGIN! INTEGER’ |; 
dl T:=T-HO;  He=H0/2; 
1FOR? zet "STEP! 1 "UNTIL? N *DO’ 
® 1BEGIN? STORELI,41:=STORELI,51; 
3 STORELI ,112=STORELI, 21 
3 "END; 
8 1 [FH 'NOT GREATER HMIN *THEN” 
2 "BEGIN? Ha=HMIN; BB3="FALSE" 
a, END’; 
ei) PRINT(* "VERKLEIN H WORDT *?,H,K); 
vs) DDz=AÂz="FALSE"; "GO TO’ P 
5 1END” 
5 "END; AAs=(KRIT ?LESS? MIN) ?AND? AA; 
2 


begin 


waarden 


rEND’; 


BEGIN? Ae ; 
805 VASKÓ(2,2,1oo=(1+X/(1-EXP(-T)))); 
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IREAL? X,PI; 

INTEGER? M; 

Pls=3,1415926536; 

IFOR Ne=t "STEP? 1 ?UNTIL! 6 'DO” 

STOREL1,N1:=0; 

Ha=HOz=HMINe==2/8; 

XseT:=0; 

MAX1=5i0=2/8; MINe=5i0=3/8; 

Mizo; Na=l; 

AAz=BBz= 'FALSE’; DDs= "TRUE!; 

PRINT(*"EERSTE ORDE SYSTEEM’ "); 

PRINT (*”KRITERIUM 1» ,MAX); 

INTCA+X,X, 1,1, R) 5 

ORG; MraM+l; 

Oil ER LGESS HATEN? 
T)); VASKO(2,2,T) 


VASKO 
1” T ’LESS’ 2o 'THEN’ ’GO TO’ R 


NN } 


jn EO} 


NO | | 


, 
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BIJLAGE 2. 


BEREKENING VOOR Bun 


De gebruikte integratieformule's zijn: 


A Ongekorrigeerd 


Ei n° 
Tar otk In 


Inr en S'n 


met DRA 


B Gekorrigeerd 


' Eve 
JT nT ntt 


Jaar In tk 


met Pi In 


Dl 
In 


T' nel 


Gezocht wordt naar oplossingen van de vorm: 


y' Pe en Yn= ee. Dan geldt: 


n+l_ 


«pe zr hee 
er“ ge” hap” 


dus P=l+in 


pe=l-in 


n+l 


ae 


h 


cn: 


Uit de beginvoorwaarden y'(O)=0 en y(O)=1l volgt 


zodat 
oek ne tf rn 
A zr l+ih) zie. ih) 


JE JL l+in)?+ HC L-ih) 5 


=P toog en I= Loo: 


Omdat h klein is kan gesteld worden : 


ltih = rei 


met r‚= Vrsn? 


lin = rei 


A 
met tg P‚=h 


ENERATOR, 


Ì voor "tijd 
verschuiving" 


J 


‚A 


Km me mm om ml oe ed md ml em al bd A LS U) 
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zodat nu gevonden is: 
Di Ds hen ee B 
g'=r, Sin np, J'=rg Sin np; 


rf OEE Ee 
Y =Fj COS np, y ; Cos n= sry Sin nf 


nu is h/2 klein dus lecos h/2 
h/?msin h/2 


mn 
us 
d SRD, { 


va 7 cos(ag;+h/e) 


Jorden nu de afw 


ingen van de uitkomsten voor y 


ewenste uitkomst y'=sin nh dan 


n 
5 8 H Ee 
zijn |= n_Mn? js moet zijn Ll 


B 


plitude in geval B neemt minder sterk toe. 


5 
ts =h moet zijn =h tg zh+t ig 
ze A 5 ze 


A 5 
ft En + eeeh + En 
moet zijn: Po=h 


5 


afwijkins P-be- F afwijkine Pg-hal- 


Dus in geval B zijn zowel amplitude- als fase- afwijking 
kleiner, (fas 


e afwijking be-fxfase afwijking A) 
Verder is in te zien dat in geval A de fouten in y en y' 
hetzelfde zijn, terwijl in geval B y over het gehele 


gebied jh in fase vóór loopt ov de fase van y'. 


=eÂl, 
Eeedoodadad 


BEGIN? 


3: 


5t 
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1BEGIN? 


BIJLAGE 5. Ë 
 wennennnnmenwnennnnnnn 


IREAL® T‚T2,H,HMIN, MIN, MAX; 
PINTEGERS NNN; REÄD(N); 


PARRAY” STORELI: N‚1s ci UITE1 EN]; 

BOOLEAN? AA,DD; 

”PROCEDURE” INTEin, RC,‚K,LAB); “COMMENT? DISAR; 
?VALUE? IN; rREAL? ÍN 7RC; 

1 INTEGER” Ks LABEL” LÂB; 


BEGIN? SWITCH” STAPz=1,2,3,4; ’GO TO” STAPCNNJ; 


1IF? DD *THEN* STOREEK ,1ds=IN; 
STORELK, 511=STORELK, 21-STORELK, 1 1=RC=H/2; 
*GO TO’ 5; 


„STORELK, 3ls=IN; 
STORELK, 51:=STÓREIK, 21-(STORELK,11+IN)=RC=H/4; 
GO: T07 55 


STORELK, 4 s=IN; 
STORELK 51 1=STÔRELK, 21+RC=H=(STORELK, 31-2=1N); 
’G0 T0*'5; 


PBEGIN?® ’REAL?’ KRIT,HI; HI g=RCxH; 


KRIT #=ABS (HT =H ( STORELK ,11-2=STORELK, &1+IN)); 
IF? KRIT ’GREATER? MAX 
AND’ H_*GREATER? HMIN eid 

BEGIN? AAz=DDs="FALSE'; NNe=1; He=d/2; Te=T2; 
1IF* H ?LESS* HMIN *THEN* He=HUlNs 
PRINT(**VERKLEIN H WORDT ”’,H,T2); 
TYPE(K); ’GO TO’ LAB 

END’; 
AAz=KRIT "LESS® MIN ’AND'AA; 
STORELK,512=STORELK, 21-(STORELK,11+4=STOREIK, 41 
HIN) /65 


PEND”; 


1[F? N ?EQUAL? K "THEN? 


*BEGIN? ’FOR? NNee=t ‘STEP? 1 ?UNTIL? N ’DO? 


UITENNI s=STOREENN, 23 8=STOREENN, 51; 
Taz=T8=T2+H; NNzel; 
PIF? ?NOT' DD 'THEN® DDz="TRUE? 
VELSEN PERZTAA ATHENS 
'BEGIN® He=a2xH; AAs="FALSE"; 
PRINT(*?H WORDT *?,H);” VASKO(4, 3,T) 
TEND’ ELSE’ AAg="TRUE'; 


vEND’; ’GO TO’ 6; 


OEEN Re K_ *THEN? 


BEGIN? ’INTEGER? 


END 


END” 


IFOR” Iz=t PSTEP! 1 *UNTIL? N DO’ 
UITCI1s=STORELI, 51; 

?IF? NN "EQUAL? 2 *THEN” Ta=T2+H/2; 
NNa=NN+1; 7GO TO’LAB 


DI SAR; 


READ(HMIN,MAX,UIT); 

FOR’ NNz=1 ’STEP? 1 ’UNTIL? N °DO? 
STOREENN, 23 z=U IT ENNJ; 

* |F? HMIN ?LESS’ „=6’ "THEN? 


*BEGIN® _ Tz=SQRT(SQRT(MAXY); Hrzos NNz=05 


1END’ 


FOR? Ha=ENTIER(T=2"POWER*NN) ’WHILE? H "NOT GREATER" 
B DO’ NNz=NNt1; H=H=e2?POWER* (NN); 
PELSE? Ha=HMIN; „PRÍnTC*STARTWAARDE VAN H *7,H); 
MINz=MAX/403 Tz=T22=0; NN t=1; 
AAz="FALSE"; DDz="TRUÊ"; 
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IREAL® A,TN, TwEePIEs 1 INTEGER? N5 

TWEEPIE:=6, banden Ne=o; TNt=1oe125=TWEEPIE; 

ed LR 

INT UITEEN 2 ‚P 

MEN ABS(T-ÍNS TESS H_ *THEN’ 

—uitvoer- IBEGIN® TNs=TN+to=TWEEPIE; Az= ARCTAN(-UITE21/UITC11); NeeNt10; 
E PRINT(N, 1 PERIODEN: AMPLITUDE ** „SQRT(UITC1 J=UIT E14 

wpdrachten UITC2I=OITL 21) 1, »1FASEVERSCHUIVÍNG 15, A-TAN=TWEEPIE); 

OIETENDALES SE Toor KIKENSSGORTOIP, 
F7 LENDE — ZEESER GOS TOL 


or 
hl SEND” 
n vEND” 
we 

| 

5 

a 

5 

En, 

K 

kj 

5 

, 

1 

en 


1 
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L 


E ORDE MIN 


ll 
} (arie dimensionaal) 


} IBEGIN? ”PROCEDURE' SE6; ’CODE'; 
YREAL? A,TN, ouD, TAU; 
INTEGER? q; 

> ARRAY? curd, hf AEL 
” SWITCH’ SCHAK tel, 
PROCEDURE’ voERUITS 


PBEGIN® WRITE(*? 
1); VASKO( 3,2 „UITC11); 
HE kr En 2 „ÛITE31); 
VASKO ‚UITL61 VASKO(2, 3 71); 
la=l+1; YEIJs=UITL61; TXL zet; 
1IF” | "EQUAL? 6oo 1THEN’ *GO TÔ' EIND 
C1 


nn 7 PRINT(*? Ca Ca Y TYD'"); 
8 TNt=o; TAUrzoe5j Az=l; OUDz=UITC 61; 
E, 
5 


: 
3,2 
3,1 
s ’y 
El 


variabelen 


Cl1ls=Cl2ls=Ci3ls=16; "are; 

VEORS [emt /STEPAeL SUNT IL? 6 ’DO’ 
1BEGIN® TX[lJs=oj YCIJe=2o=( 1-1) 
VEND’; liemós 


me ed 
| begin- 


en 
SP: INT(CL13,TAU, 1 
3 INT(CL21,TAU, 2 
5 INT(CL31,TAU, 3 
5 
kei 
8 
sl 
dÁ 


1 
+UITL4],10/3, 
T(UITL4I+UITE 5] 
(UITL51+UI TC 6) 5 
T+H/2-ENTIER 4 ’LESS’ H 'THEN’ VOERUIT; 
UITC61-OUD ’GREAT THEN” 
» UITL6I 'GREATER’ 30 ‘THEN’ 
s=o.l "ELSE TAUzz=o.0o2; 
PIF? T 'GREATER® TN 'THEN* 'GO TO’ SCHAKLQI 
IEND’ 'ELSE’ '|F? UITC61- ‘GREATER’ 30 ’THEN? 
TAU:=zo.5 'ELSE’ TAU:=o.02; 'GO TO’ 4; 


AB, 
ze); 
‚P 
INT(=UIT 1 J=UI TL 
/3 
[5 
6 


or 


5 
: 
2 
E 


C[QIs=-CLQI; Q2=Q+1; TN:=T+A; 
VOERUIT; ’GO TO’ 4; 


C[3ls==-C[3]; Qe=i; TNs=T+A; 
VOERUIT; 


OUDz=UIT[61; "IF! T 'LESS” 6o *THEN’ ’GO TO’ P; 
SE6(I,1,Y,TX, TRUE”); 


Gl 


BIJLAGE 6. SIMULs 


KE PRINT(**KRITERLUM?*,MAX); 
IBEGIN® ?REAL? Di,D2,NUo,NU1, NU2, NU3,W1,W2,VX1,VX2,VY1,VY2, 
en PSIA,PSIE,VA,VE, VS,F6,Fo;A; 
PARRAY?® P‚QLo:231,CÍ1:271, AVE, AVA, AT[1:415 
En PREAL? PROCEDURE” NULPUNT(X,Y); 
PARRAY? XY; 
"BEGIN® 'REAL* A,OPL; 


Í 
ef PARRAY? B[L1:41; 

ie Ee 1 INTEGER? M‚N; 

zi ig OPL 1=0; 
85 1FOR? Ne=1 "STEP’ 1 *UNTIL® 4 DO’ 

5 EE: IBEGIN® 'FOR* Mz=2,3,4 'DO’ 

| s 6 PBEGIN® BL1lr=XEND; 
EK: EN BEMJz= "IF” M "EQUAL? N ’THEN? X[11 'ELSE’ XIMJ; 
1END’; 

d 2 At=YEND; 

B. el OPL s=OPL-BL21=BL31=BL4I=A/((BL11-BL2I)= 
d (BL11-BL31)=(BL13-BL4I)); 

— je, PEND’; 

E El NULPUNT #=0PL; 
> "END? PROCEDURE NULPUNT; 

5 


PREAL” ’PROCEDURE’ F12(X); 
— 'VALUE' X; 'REAL” X; 
YBEGIN® REAL’ Zy; 

Zr=ABS(X); 


, IF! Z ?NOT GREATER® C[1] THEN’ Ys=Z=<C[21/C[31 'ELSE' 
IF” Z NOT GREATER? 4 'THEN? Ya=Cl4l=Z=Z-C[5]=Z+ 
a C[6] ’ELSE’ 


’IF” Z ’NOT GREATER’ C[7] 'THEN’ Ya=C[8)=Z-C[9] ELSE’ 
LF” Z *NOT GREATER’ C[10]l THEN? Ys=CCtil=Z=Cl121 
A ke SELSE 


pen 


8 1[F* Z 'NOT GREATER” C[131 'THEN’ Yz=C{141=Z2-C[151 
E. re ELSE’ Y3=5; 
® F122=SlGN(XS=Y 
5 D ZEND’; 
Ë. D 1REAL? ’PROCEDURE' Fo6(X); 
© IREAL’ X; 
= B 'BEGIN’ 
3 
5 


1IF! Z 'NÔT GREATER' C[16] 'THEN’ Ya=Z/CL171 "ELSE 
ae )IF? Z "NOT GREATER” C[18J ’THEN’ Ys=Cl19}=ZeZ- 
GEZol=ZtCL2 LI ZELSER 


— El IF! Z *NOT GREATER? C[22] *THEN? Ys=C[231=Z-C[ 24) 
Sk "ELSE? 

E rd 1IF! Z *NOT GREATER’ C[251 "THEN' Ya=CL261=Z- 

4 Ei C[271 'ELSE* Ys=0.7; 
Ë Foós=SlGN(X)=Y=105 


(vervolg 


inlezen 


potentiometer 


[é) 


andere 
beginvoorwaarden 


BEGIN’ 
END’ 
_ BEGIN? 
relais 
met ZEND’; 
if BEGIN? 
ein 
Peren END’ 
Eras 'BEGIN 
EPEN 


integratoren 


ZeOeZe 


“Di t=D23=0; 
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READ(C,PC11,QC191,QL213,QL41,PC9I,PC231,QL12), 
P[21],P[20],VS,P[41,P[2],QL61); 

Plols=P[1]; Q[ 20 P[20]; 

QLols=QCBIE=QCI [18):=Q[191; 

P[3]s=Ql211; QC31s=QL41; 

P[L8)s=P[9 QL1ols=QL121; 
P[tols=P[1213=Q[22l:=P[23); 


VROR? Aumtd STEP? 4 ZUNT IE? ASZ DO 
ATCA]s=AVALAI 3=AVELA] z=A; 

PS Azs=to=UITL11; PSIEs=UIT[ 2]; 
VAs=UIT[ 3); VEs=UIT[4); 

P[15]:= VE-to=PSIE=P[23]-D2; 

QL15ls= VErLo=“PSIE=Q[L221-D1; 
NUor=-to=Pltol=F12(P[15]) — PL4]=PSIA; 
NUts=to=QL61=VA+PL21=PSI A-P[21]=PSlExlo; 
NU2s=1o=Q[6}=VA-P[2]=PSIA-P[21]=PSI Esto; 
NU 3s=-P[4I=PSI Amto=P[12}=F12(Q[151); 
ERE 
W2s=ABS(NUo)+ABS(NU2); 

|F? Wi ’LESS’ to =P[9] ’THEN? 
VXis=NU3; VY1s=NU1 

VELSEN 

VX1 3=-QL9 1=NU 3/W1; 
VY12=-Q[81=NU1 /W1 


wiee Wa ZEESS LOP CB ATHENS 
VXaz=NUo; VY2s=NU2 

HEESEN 

VX2e=QL191=NUo/W2; 
VY2r=-Q[181=NU2/w2 


Ee 
Foz=Fo6(P[151); 
INT(P[2ol=Fo+Qitol=VX2+QL12J=VX1-QL201=F6-QL211=VY1+ 
PL3]=VY2,=to,1,START); 

INT(PSIA, =1o,2,START); 

et nn de VEL 
INTÚ=1o=VA, 10,4,START); 


_ 


mn mn a 
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vervo ly 


SFOR’ As=3,2,1 ’DO' 
ATLA+11 
AVELA+1 


AVALA+1]s=AVACAI; 
,END'; : 
ATL1]s=T; 
AVEL1Js=UITL4); 
AVAL1]s=UITC31; 
YIF’ T_ ’GREATER® oat *THEN? 
B BEGIN? *IF* SIGN(AVEL21) ’NOT EQUAL? SIGN(AVEL31) 'THEN 
uitvoer nulpunt PRINT(*?NULPUNT T="?,NULPUNT(AVE,AT)); 
uitvoer i 1|F? SIGN(AVAL21) 'NÔT EQUAL’ SIÓN(AVÁC31) *THEN? 
of PRINT(? MAXIMUM T="" 


BNET „ AT), NULPUNT (AVA, AVE); 
TEND’;* IF*T+H/D-ENTLER( 2ont+5eli)/20 'LESS’'H "THEN? 
BEGIN? VASKO(1,3,T); VASKO(1,4,UITL4I) 

1END'; 'IF? T 'LÊS$? 1,5 'THENS 9GO TO’ START 


t.b.v. DISAR: 
BODO 00) 29 

tebev. functie neratoren: 
„45,.15,2.45,eÂl, 017, 
12,96,4,26,3.8461,13,58461,4,37, 
6,3636,24,30909,4.5, 11. 53846, 46.92307, 
NE 25,…21534, 1e 5227, 
2e 78136, eh 56,1,2903,5.2838,4.7, 
„71428, „65714, 

votenti enen 


„5123,.2,=5142,.8012,19+2 
dien, 
25,.315,625,1 


